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REMARKS 

1.0  SUMMARY 


Large-disturbance  equations  of  airplane  notion  suitable  for  use 
with  modern  digital  or  analog  computers  are  derived  from  the 
general  equations  of  rigid-bod/  notion.  The  procedure  follcued 
is  to  find  solutions  in  each  case  for  the  steady^fllght  conditions 
about  which  the  desired  notions  .eptesent  perturbations.  The 
variablec  of  motion  are  expended  to  Include  these  steady-flight 
values,  and  the  steady-flight  solutions  are  applied,  in  order  to 
eliminate  the  steady-flight  forces  and  moments. 

\ 

Two  specific  cases  are  carried  out  to  final  forn.  These  are: 

(1)  Six  degree -of- freedom  equations,  based  on  arbitrary 
body  axes,  with  small  longitudinal  velocity,  sideslip, 
and  angle  of  attack  perturbations. 

(2)  Five  degree-of-freedom  (constant  longitudinal  velocity' 
equations,  based  on  principal  axes,  with  small  sideslip 
and  angle  of  attack  perturbations,  and  with  the  effects 
of  gravity  simplified. 

The  first  set  of  equations  is  shown  to  be  especially  suitable 
for  use  in  fire-control  or  tracking  studies,  while  the  second 
set  of  equations  is  intended  for  use  in  calculations  involving 
r-pid  rolling  maneuvers. 

Auxiliary  equations  are  developed  for  the  conversion  of  wind- 
tunnel  data  to  body-axis  stability  derivatives  and  for  the 
iterations  required  to  obtain  numerical  solutions  of  the  steady - 
fjight  equations  of  motion . Also,  as  an  aid  in  calculating  and 
interpreting  instrument  readings  in  conjunction  with  large-dis- 
turbance motions,  the  readings  of  attitude-,  velocity-,  ar.d 
acceleration-measuring  instruments  are  derived. 


DOUGLAS  A»CtA>l  COMPANY,  INC.  HVGUNDO  DIVISION  tl  SCGUNDO.  0*Uf0«NIA 


Page  2 
ES  17935 


2.0  TABLE  OF  COHTSKTS 

Paso 


1.0 

Sunnary 

1 

2.0 

Table  of  Contents 

2 

3.° 

Introduction 

4 

4.0 

Reference  Amec,  Symbols,  and  Definitions 

5 

4.1 

Reference  Axes 

5 

4.2 

Symbols 

5 

4.3 

Definitions 

10 

5.0 

Development  of  the  Equations  of  Motion 

11 

5.1 

General  Equations  of  Airplane  Motion 

11 

5.2 

Steady-Flight  Equations  of  Motion 

• 13 

5.2.1 

Steady  Pitching  Symmetric  Flight 

13 

5.2.2 

Steady'  Straight  Symmetric  Flight 

14 

5-3 

Expansions  of  the  Variables  to  Include 

Initial  Conditions 

14 

6.0 

Special  Forms  of  the  Equations  of  Motion 

19 

6.1 

Equations  Suitable  for  Fire-Control  Studies 

19 

6.2 

Equations  Suitable  for  Calculation  of  Rapid 

Rolling  Maneuvers 

20 

7-0 

Stability  Derivatives  on  Body  Axes 

21 

7.1 

Longitudinal  Velocity  Derivatives 

21 

7.2 

Angle  of  Attach  Derivatives 

24 

7-3 

Sideslip  Derivatives 

• 24 

7.4 

Symmetric  Rotary  Derivatives 

25 

7o 

Asymmetric  Rotary  Derivatives 

25 

7.6 

Control  Surface  Derivatives 

25 

8.0 

Development  of  Auxiliary  Equations 

26 

8.1 

Initial  Condition.  Iterations 

26 

G.2 

Readings  of  Attitude -Measuring  Instruments 

26 

n 
**  • 

x.cnii.ju*0r>  >.>j.  v co-Ou jl instruments 

27 

8.4 

Readings  of  Aece] e rat  ion-I leas . iring  Instruments 

28 

8.4.1 

Linear-Typo  Accelerometer 

31 

8.4.2 

Pendulum -Type  Accolerome ter 

32 

8.4.3 

Pendulum -Tyne  Angular  Accelerometer 

34 

9.0 

References 

35 

o 


h|»  3 

K 17935 


BOjlftM  CQM»MW,  INC  ft  MOUNPO  0»V«K»N  H »KX*«0.  CM*0*MA 


TABLES 

£US 

I 

Equation*  of  Airplane  Motion  Suitable  for 
Fire-Control  Studies 

A-2 

11 

Equations  ol  Airpluna  Jiotion  Suitable  for 
Calculation  of  Rapid  Rolling  Maneuvers 

A-3 

III 

Longitudinal  Velocity  Derivatives  on  Body  Axes 

A4 

IV 

Angle  of  Attack  Derivatives  on  Bcdy  Axes 

A-5 

V 

Sideslip  Derivatives  on  Body  Axes 

A-6 

VI 

Symmetric  RotBr;  Derivatives  on  Body  Axes 

A-7 

YII 

Asymmetric  Rotary  Derivatives  on  Body  Axes 

A«£ 

VIII 

Control  Surface  Derivatives  on  Body  Axes 

A-9 

FIGURES 

1. 

Airplane  Body  Axes 

A- 10 

2. 

Angular  Relationships  in  Plane  of  Symmetry 

A-U 

^ • 

angles  of  Attack  and  Sideslip 

A- 12 

4. 

Orientation  Angles  of  Body  Axes 

A-13 

5- 

Aerodynamic  Force  Resolutions 

A*  lb 

6. 

General  Instrument  Axes 

A- 15 

n * 

Accelerometer  Orientations  Relative  to 

Instrument  Axes 

A-16 

APPENDIX 

A. 

The  Reading  of  a Linear  Acceleroscter  at 
the  Airplane's  Center  of  Gravity 

A-l 

co— we.  mc  n hounoo  omson  « mO'Moo.  coupo— 


3.0  HTTRODUCTIOW 


The  equatiors  of  rigid  body  notion  referred  to  body  axes  de- 
veloped by  the  Swiss  mathematician  Leonhard  Euler  are  the  basis 
for  all  airplane  dynaalc  stability  theory.  These  equations  have 
been  most  widely  used  In  aeronautics  in  their  linearized  form, 
and  concepts  such  as  the  period  and  time  to  damp  to  half-aaplltude 
of  an  oscillation,  the  neutral  point,  and  the  transfer  function 
are  rigorously  defined  only  in  terms  of  the  linearized  equations. 

Linearization  of  Euler's  equations,  as  proposed  by  E.  J.  Bouth, 
requires  that  the  motions  described  be  restricted  to  3mall  per- 
turbations about  an  initial  condition  of  steady  notion.  An 
important  slssplifi  cation  in  the  aeronautical  field  occurs  when 
the  initial  steady  motion  is  a case  of  "symtetric  ' flight,  in 
which  the  airplane's  pi  one  of  symmetry  remains  fixed  in  a ver- 
tical positio;..  Small  disturbances  from  steady  symmetric  flight 
are  described  by  two  independent  sets  of  three  simultaneous 
equations  of  notion,  instead  of  by  the  general  single  set  of 
six  sijaultoneous  equations.  In  the  general  llneerized  case, 
however,  the  initial  steady  motion  about  which  small  perturba- 
tions occur  can  be  a severe  maneuver,  such  a?  a steeply  tanked 
turn,  a rapid  roll,  cr  a spiral  dive.  It  is  only  in  cases 
where  the  perturbations  or  disturbances  themselves  must  be  allowed 
to  be  large  that  the  linearized  equal  1oji3  of  notion  are  inappli- 
cable. 

The  development  of  certain  special  non-linear  forms  of  the 
equations  of  airplane  notion  which  apply  to  large  disturbances 
from  steady  flight  is  the  subject  of  this  report.  Specific 
applications  vnich  have  been  considered  are: 

(1)  Fire-control  studies,  where  the  airplane  soy  be  re- 
quired to  perform  extreme  maneuvers  while  tracking  a 
target 

(2)  Caic-iiTcion  A'  r^vfu  dTxfng'm'ir. silvers  and  similar 
flight  conditions  possible  with  modern  airplanes 

The  inadequacy  of  the  linearized  c rr-ull-pcrturbn tl on  equations 
of  motion  for  either  of  the  e applications  has  been  established. 

In  the  expression  of  the  equations  of  this  report  ir  fern  suita- 
ble for  machine  computation,  iongitudinul  velocity  and  angle 
of  sideslip  and  attack  disturbances  are  restricted  to  small 
values,  as  in  the  linearized  theory. 
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4.0  REFERENCE  AXES  t SYMBOLS  AND  DEFINITIONS 


4.1  Reference  Axe; 

The  large-disturbance  equations  of  airplane  notion  developed  in 
this  report  are  baaed  on  systems  of  axes  which  are  fixed  in  the 
airplnne,  or  body  axes.  It  is  often  expedient  to  consider  a 
special  case  of  body  axes  called  principal  axes,  in  which  the 
X-,  Y-,  end  Z-*ax'es  coincide  with  principal  exec  of  inertle. 

For  any  set  of  body  axes,  however,  the  following  definitions 
are  used.  The.  oxes  are  orthogonal,  having  their  origins  at 
the  airplane's  center  of  gravity.  The  X-axis  iie6  within  the 
airplane's  plane  of  symmetry,  and  is  positive  forward.  The  Z- 
axi3  is  also  in  the  plane  of  symmetry,  and  is  positive  towards 
the  bottom  of  the  airplane.  The  1-axis  is  perpendicular  to 
both  X and  Z,  and  is  positive  towards  the  right  wing.  Positive 
senses  of  quantities  referred  to  body  axes  are  governed  by  the 
positive  directions  of  the  axes  themselves,  following  the  right- 
hand  rule  in  the  case  of  angular  velocities  and  moments.  Body 
axe3  ere  illustrated  in  Fi^re  1 . 

Basic  aerodynamic  data  used  in  the  equations  of  this  report 
are  assumed  to  have  been  obtained  on  wind-tunnel  stability 
axes.  The  wind-tunnel  stability  axis  system  is  defined  identi- 
cally to  airplane  body  axes  except  that  the  Z-3tabtlity  axis 
is  perpendicular  to  the  relative  wind,  or  approximately  to  the 
wind-tunnel  axis.  In  addition,  the  origin  of  wind-tunnel  sta- 
bility axes  may  not  coincide  with  the  assumed  airplane  center 
of  gravity. 

4.2  Symbols 

Length  is  measured  in  feet,  mass  in  slugs,  time  in  seconds, 
and  angles  in  radiens,  unless  specified  otherwise. 

b * Wing  span 

c - Vi iii*  ciis-u  Merouynamic  Chord} 

Cg,  C^;  3 Drag,  lift,  and  thrust  coefficients.  Cp  =>  D/Sq 

Cl  * L/Sq 
Tg  « T/Sq 


CIv,  Cr,  Cr  Cj  , Ct  • Moment  of  Inertia  and  product  of 
* nf  2 X2  inertia  coefficient*. 

cix  • Vs»ib 
CIy  * W S(liU 
ciy'  - ig/tsq^b 

CIZ  * W 
cixz  ■ J0t2/S,ib 

Cjj  CB,  Cn  * Rolling,  pitching,  and  yawing  moment  coefficients. 

C|  3 L/5(j' » 

CB  - M/Sqc 
Cn  - N/Sqb 

Cy,  Cy,  Cz  * Longitudinal,  side,  and  normal  foxce  coefficients. 

Cx  - X/Sq 
Cy  • Y/Sq 
CZ  - Z/Sq 

C(  >u- ' C(  V )a  * >/0U ’ ’ *C(  d°(  )/9j 

C(  )i'  C(  )y  C(  )•  * *>C(  y^('>  3C(  )/3r 

c(  )p,  C(  } . C »0C(  j/cHpb^),  dC(  )/0(qc/2V1), 

dC(  j/dtrt/^) 


C/  ^ ? C.  . , etc. 

I '3a  ^7*  V '« 

d «*  Aeceleroxeter  dumping  constant 
D,  L,  T ■ Drug,  lift,  and  n=:t  thrust 
f ■ Accelerometer  pendulum  length 
[g]  * Gyro  matrix 
g ■ Acceleration  of  gravity 


H = Endin'.  uiii^ujHr  mor.v.iitum  (defined  us  positive  for  right-hand 
rotation,  viewed  from  the  rear) 


•Sq  - *T/«V 

u,  r,  w » Disturbance  longitudinal,  lateral,  and  normal  linear 
velocities  of  body  axes 

U,  V,  W ■ Longitudinal,  lateral,  and  normal  linear  velocities 
of  body  axes 

u'  ■ Dimensionless  small-disturbance  velocity  variable,  u/V^ 


V • Oross  weight  of  airplane,  mg 

x,  y,  a - Distances  along  X,  Y and  *'•  body  (airplane)  axes 

X,  Y,  2 ■ Longitudinal,  side,  and  normal  aerodynamic  forces 
along  body  axes 

X(  )'  y(  )’  r(  ) “ B<*^  axes  stafcility  derivatives 


x/  \ " 3 X 1 
1 } 


y(  ) • -14  JL 

1 dTT  qxS 


5T7  qxS 


Perpendicular  distance  of  thrust  axis  below  center  of 
gravity 
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i,  J,  k - Unit  vectors  along  X,  Y and  Z body  (airplane)  axes 

*X'  JY'  I2»  ^ ■ Moments  of  inertia  about  the  X-  Y-  and  2-  axes, 
and  product  of  inertia  about  X and  2. 

I - Accelerometer  pendulum  polar  moment  of  inertia 

if  • Angle  of  incidence  of  thrust  axis  to  X-  axis. 

k - Accelerometer  spring  constant 

*(  )»“()»  n{  ) ■ ®04y  axes  stability  derivatives 

1/  > - OL  i 

( J an  q^b 


■/  \ - Dm  i 

( 3 dT7 


n(  ) - «>N  1 

d(  ) q^Sb 

1,  n,  n n Unit  vectors  along  *I#  Yj,  and  2j  instrument  exes 
L,  M,  N ■ Rolling,  pitching,  and  paving  aerodynamic  moments 

[L]  « Orientation  matrix 
n ■ Load  factor 
n « Mass  of  airplane 

* Accelerometer  celomic  mass 

“ ^ --  rPrVnt^^rP,.--:.  r.7,]  _.id  javiiig  aaigulttx- 

velocities 

P'  Q,  P.  = Rolling,  pitching,  and  yavlng  angular  velocities  of 
body  axes 

q « Dynamic  pressure,  (P/P)  V“ 

r = Radius  vector  of  instrument  location  to  airplane's  center 
Of  gravity 

s “ Op-***:  tor  d/dt 

t,  u = '-'nit  vectors  relating  to  accelerometer  (See  Figure  7) 
S • Wing  :;r  \ 
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B' , a'  • Dimensionless  small-disturbance  velocity  variables. 

B*  ■ v/Vx 
o’  - v/Vj^ 

B,  a ■ Sideslip  angle  and  angle  of  attack. 

B - *m‘1  V/V 
a « tan**  V/U 
- Plight  path  angle 
& « Control  surface  disturbance  angle 

€ - Angle  of  nose-up  incidence  of  fuselage  reference  line  relative 
to  X - principal  axis 

J » Accelerometer  damping  ratio 

^ - Angle  of  attack  of  X - principal  axis 

^ - Air  density,  also  tilt  angle  of  sensitive  direction  of 
accelerometer  to  the  horizontal 

T • Time  factor,  m/^^j 

i,  9,  • « Orientation  angles  of  reference  axes  to  arbitrary  earth 

% axes 

^ » Disturbance  values  of  the  orientation  angles 

Jl  r Total  angular  velocity 

<t>n  « '.ocelerometer  undamped  natural  frequency 
Subscripts 

a,  e,  r Aileron,  elevator,  rudder 
s = Accelerometer  sensitive  direction 

t * Normal  to  pendulum  accelerometer  support  axis  and  rensitive 
directions 

u » Pendulum  accelerometer  support  axis  direction 
I » Instrument  or  instrument  axes 

0 » applied  force  or  moment 

1 >■  Gtrad.v  flight 
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*»  .3  Peflnltlor-* 

Steady  p.lght  1b  defined  as  flight  with  zero  rates  of  change 
or  the  linear  and  angular  velocity  variables,  or  with 

0-V»W«P*Q*R-O.  Steady  sideslips  and  steady  level  turns  or 
steady  helical  turns  about  a vertical  space  axis  are  possible 
steady  flight  conditions.  Steady  pitching  flight  lo  more 
properly  referred  to  as  ** quasi-steady" , since  0 or  V>  cannot 
both  retain  zero  over  a long  period  of  tine  If  QHO. 

Straight  Flight  is  defined  as  flight  with  zero  values  of  the 
angular  velocity  variables  P,  Q,  and  R.  Steady  sideslips 
and  dives  or  climbs  with  or  without  longitudinal  acceleration 
are  sone  straight  flight  conditions. 

Symmetric  Flight  Is  defined  as  flight  with  a fixed,  vertical, 
position  in  space  of  the  plane  of  aynnetry.  In  symmetric 
flight,  the  variables  P,  R,  *,  and  V remain  zero.  Wings* 
laterally* level  dives,  climbs,  pullups  and  puehdovns,  with 
zero  sideslip,  are  »ome  symmetric  flight  conditions. 

Asymmetric  Plight  Is  defined  as  flight  in  which  the  plane  of 
symmetry  does  not  remain  in  a fixed  vertical  position.  Some 
or  all  of  the  variables  of  motion  P,  P.,  ♦,  or  V may  be  expected 
to  be  other  than  zero  in  asytnitric  flight.  Examples  of 
asymmetric  flight  include  sideslips,  rolls,  and  turns. 


row  }8-2*fl* 
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>.0  DtoVKIOFMElTT  OF  THE  EQUATIONS  OF  MOTION 


5 . 1 Oengral  Equations  of  Airplane  Motion 


The  general  equations  of  airplane  notion  are  a set  of  nine 
simultaneous  non-linear  differential  equations,  listed  here 
as  Equations  (l)  to  (3): 


Applied 

Aerodynamic  Gravitational 
Poreas  Forces 

X - ng  sin  8 

Y + ng  cos  8 sin  ♦ 


♦ ng  cos  8 cos  * » ■ 


Rates  of  Change  of 
Linear 
Momentum 


(U-VR+VQ) 

(v-wtmjf) 

(W-OQ+VP) 


Applied 

Aerodynamic 

Moments 


Rates  of  Change  of 
Angular  Momentum, 
Engine  Stepped 


Rates  of  Change 
of  Engine  Angular 
Momentum 


- IjfctR+PQMVIy)®  QHe  sin  I?  ' 

M ♦ IxsO,2-',2Hlx-^)pR  ♦ slniT+  R^  cosiA^) 

" CRe  cob  i^.  I 


Rates  of  Change 
of  Orientation 
Angles 


Function#  of 
Angular  Velocities 


y * (Q  sin  * + R cos  «)  ton  0 


■ Q cos  • - R aln  ♦ 


Equations  (2)  .pply  to  airplane*  in vlng  mirror  lyrmetry  about  the 
X-Z  plane  except  for  engine  angular  momentum.  Consequently, 
the  effects  of  rudder  and  aileron  deflections  cn  the  inertia 
parameters  are  neglected.  The  aerodynamic  fore-?#  and  moments 
in  Equations  and  (2)  are  usually  found  to  depend  upon  the 
lli«r  %r  and  anguj-A*/  velocities  and  their  derivatives,  the  control 
surface  angles  and  chelr  derivatives,  air  density,  Mach  number, 
and  possibly  other  factors  as  veil.  These  functional  dependences 
are  represented  by  Equation  (4). 


ro«M  jo-jsor 
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5.2  8teady-night  Equations  of  lotion 

fbe  steady-  flight  out  of  the  general  equations  of  airplane 
■otlop  la  uaeful  for  establishing  Initial  conditions  for  the 
numerical  aolution  of  ttaa  general  equations,  as  veil  as  for 
ealoulstions  involving  steady  flight  only.  In  the  aost  general 
expression  tbs  steady- flight  equations  of  airplane  notion  a*-# 
obtained  froo  Equations  (l)  and  (2)  by  eliminating  tbs  terns 
involving  derivatives  of  the  v&rlableeOffVjVjP,^,?),  and 
eeploying  the  steady-flight  subscript  1 throughout. 


Two  special  cases  of  the  steady*  flight  equations  of  airplane 
notion  are  of  particular  interest,  and  these  cases  are 
diacussed  in  the  following  sections. 


5*2.1  Steady  Pitching  Symmetric  Flight 


Boiling-pullout  or  pushdown  saneuvers  any  be  performed  by 
applying  abrupt  aileron  deflections  while  in  steady  symetrlc 
pull up  or  pushdowns . Using  the  dimensionless  forms  of  Equations 
(l)  and  (2),  with  o obtained  from  Equation  (5),  the  steady- flight 
equations  which  apply  are: 


Cxi  ein  (*x  ♦ eq)  » 2TQi  sin  aq 

Cy  • 0 
U T 

009  (*1  ♦ • -STQj,  cos  (X^ 


7) 


y 


The  steady  pitching  velocity  is  often  expressed  in  terns  of 
the  normal  load  factor  n^.  Load  factor  la  defined  as  the  ratio 

of  applied  aerodynamic  force  In  the  specified  direction  to  the 
gross  weight.  Thus: 

°z  “ * y (6) 

The  negative  sign  in  Equation  (8)  is  U3ed  to  obtain  agreement 
with  the  usual  sign  convention.  Substituting  Equation  (6)  into 
the  third  of  Equations  (6)  and  solving  for 

[ "Zi  ' co*  <ri  4 *j] 


" ifjcoeaq 


(9) 


fi 

i f 

* 
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The  maeileal  solutions  of  Equations  (6),  (7)  and  (9) 
discussed  <n  Section  8.1.  Because  of  the  arbitrary  functional 
dependence  of  Cg  and  Cj  on  angle  of  attack,  this  eolution 

la  found  by  iteration. 


Steady  straight  symmetric  flight  may  be  used  as  a reference 
for  studying  almost  any  maneuver,  using  the  large -disturbance 
equations  of  airplane  motion.  However,  because  of  practical 
limitations  in  the  way  of  inforumtlon  storage  on  either  digital 
or  analog  computers,  it  will  be  found  expedient  to  establish 
Initial  or  reference  conditions  which  are  as  close  to  the  desired 
maneuver  as  is  possible.  For  example,  although  steady  straight 
symmetric  flight  may  be  used  as  a reference  for  studying  the  roll* 
log  pullout  maneuver,  initial  conditions  which  include  the  steady 
pitching  portion  of  the  maneuver  (as  in  5 .2.1)  eliminate  the 
necessity  for  machine  ccaputction  of  that  much  of  the  problem. 

Where  steady  straight  symmetric  flight  is  chosen  as  the  Initial 
steady  flight  condition,  the  applicable  equations  of  motion  are 
Equations  (6)  and  (7),  with  * 0. 

Expansions  of  the  Variable n to  Include  Initial  Conditions 


Ihe  steady-flight  portions  of  the  general  equations  of  airplane 
notion  can  be  eliminated  from  the  equations.  As  previously 
indicated,  this  will  reduce  the  amount  of  information  storage 
needed  for  nunerlcal  solutions  of  the  equations.  The  steady- 
flight  portion  of  the  equations  may  be  eliminated  if  the  variables 
of  the  motion  are  expanded  to  Include  Initial  conditions.  This 
is  done  in  the  following  equations  for  the  attitude  angles,  the 
angular  velocities  and  i ccelerations,  and  the  linesr  velocities 
and  accelerations.  It  should  be  appreciated  that  in  these 
equations  the  disturbance  quantities  indicated  by  lover-ceee 
symbols  or  by  the  prefix  A are  not  necessarily  small  disturbances, 
as  they  are  in  Reference  4. 

0 ■ 0j_  ♦ & (10) 


sin  0 • sin  6X  cos  &+  cos  sin 
cos  0 •»  cos  0i  cos  <*-  r«n  0,  sin  & 


♦ • ♦ ip 


H<  NWi 


jo-*: 

■51) 
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V » u • -(VL  > 4T)  cos  A3  (sin  cos  Aa  ♦ cos  slnAa)Aa 
-(Vx  ♦ &)  sin  A3  (cos  Ox  cos  Aa-sin  elnAa)A3 
♦ df  cos  A3  (cos  cos  Aa  • sin  Oj,  sin  Aq) 

v - V - (Vx  4 /y)  cos  A3  a£  4 Af  sin  A3 

W • w ■ (Vx  ♦ 40  cos  A3  (cos  cos  Ac  - sin  c*x  sin  Aa)A& 

•(Vx  ♦ 40  sin  A0  (sin  cos  Au  4 cos  sin  Aa)A3 

4 AV  cos  ^ (sin  cs,  cos  Aa  * cos  »in  Aa) 


Equations  (18)  to  (ZO)  express  linear  velocities  and  accelerations 
along  airplane  body  axes  in  teres  of  the  large-disturbance  variables 
Ap,  A3,  and  Aa.  Aerodynamic  flow  breakdown  always  acts  to 
limit  the  practical  ranges  of  these  variables,  particularly  the 
latter  two.  Ibis  makes  it  expedient  to  use  the  small -disturbance 
linear  velocity  variables  if , 3',  and  a'  in  place  of  4f.  A3,  and 
Aa  respectively,  even  though  the  angular  displacement  and  velocity 
disturbance  variables  *>^’,4*,  p,  q,  and  r are  allowed  to  take 
on  large  values.  The  small-disturbance  linear  velocity  variables 
are  defined  as: 


Following;  2.M.  Jones  on  page  146  of  Reference  4,  the  total  velocity 
change  A|f  is  expressed  in  terms  of  u’  and  a'  as  follows: 

From  Figure  3:  * 

V2  - U2  4 V2  4 WS  (22) 

Differentiating: 

2Vx  dV  “ 2Ux  u 4 2Vx  v 4 2V}W 
Since  we  assume  that  Vx  ■»  0: 

--  u'  cos  c*x  4 a'  sin  ax  (23) 

From  Equation  (lfy) 

A3  = cin"^  v 

Vi-4’ 


F4M  J0-? 
-*-SI»  ... 
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8inc«  both  v and  d*  ere  mil  the  elne  equals  the  angle  and 
Of  !•  negligible  coopered  withVx*  Thus,  from  Equation*  (21) 

£$  m p'  (2*) 

Proa  Equations  (l6)  and  (19 )’ 
ta  ■ a - 

. w1!^  - un-1 

Ux+U  Ux 

«X+*  «x 

• tan"1  Ux  ♦ u * 

, (Wi+w)  wT 
1 +(Oi+u)  Ul 

2 2 2 

Simplifying  and  substituting ^x  **or  ^1  + V1  > 

cos  for  OxyVi  » ^ *1®  °l  for  wiATi  » 

Aa  - tan"1  v COB  " u sln  Q1 
*Vx  ♦ u cos  oq*v  sindt^ 

v and  u sre  both  assumed  smell,  and  negligible  compared  with 
Thus,  ftrom  Equations (21):  , 

Act  » o'  cos  ax  - u'  sin  ax  (25) 

Note  finally  that  since  P*  and  a'  are  small, 


sin  P'  - p' 
sin  a*  « a* 


cos  P'  ■»  1.0 
cos  a*  * 1.0 


and  products  of  u',  0' , a',  and  their  derivatives  are  negligible. 

Equations  (23)  to  (26)  nay  now  oe  substituted  into  Equations  (19) 
and  (20),  leading  to: 

U - Vx  (u'  ♦ cos  ax) 

V - Yl  p27* 

W - ^x  (a'  ♦ oin  ax)  J 

0 - “Vx  «'  "1 

V • \i  p’  p28) 

W - J 


commmt.  mc  n wouMOQ  wva<ON  atmiNOft  wawt 


yw*, 


Equations 
Equations 
(vlth  Vx  ■ 


can  of  course  be  obtained  directly  1. 
using  t^ese  relationships  froa  Figure  j 


cos 

sin  - Wly^i 


In  onler  to  obtain  numerical  solutions  of  the  equations  of 
airplane  motion  the  aerodynamic  forces  and  moments  must  be 
expressed  in  terns  of  the  linear  and  angular  velocity  variables 
of  motion.  In  the  general  case  of  large  disturbances,  these 
expressions,  symbolized  in  Equations  (4),  can  be  exceedingly 
non- linear  and  complex.  These  difficulties  arise  from  aerodynamic 
flow  separation  at  large  relative  flo*  angles  and  at  high  subsonic 
Msch  numbers,  from  tbe  Interacting  effects  of  large  values  of  the 
variables  like  sideslip  and  angle  of  attack,  and  from  possible 
aerodynamic  hysteresis  effects.  Where  aerodynamic  non-llnearltles 
can  be  disregarded,  Equations  (30)  are  Illustrative  of  the  exparslon 
of  the  aerodynamic  forces  and  moments.  These  equations  are  and* 
dimensionless  to  conform  vlth  Equations  (1A)  and  (2A),  and  the 
small-disturbance  linear  velocity  variables  are  employed. 


JL  . cx,  * -i-  (?5-  «•  . c . 

q1S  X1  qjS  Vu'  do 


-i—  < CY,  +■ 

qjS  H 


c 


I 

I 
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6.0  SPECIAL  FORMS  OF  THE  E^UmOHS  OF  HOT70B 

Vbe  ger_-ral  equations  of  airplane  notion  for  arbitrary  body  axes  are 
given  with  reasonable  generality  in  Equations  (l)  to  (5).  Several 
special  forms  of  these  equations  are  developed  in  this  section, 
making  use  of  a series  of  approximations  which  are  designed  to 
simplify  the  equations  as  far  as  possible,  consistent  with  their 
Intended  uses. 

6.1  Equations  Suitable  for  Tracking  Studies 

In  studying  the  closed-loop  dynamics  of  an  automatic  tracking 
or  fire -control  loop  it  is  typical  to  calculate  the  notions 
of  the  attacking  airplane  from  the  initiation  of  automatic 
tracking  to  t.\e  firing  point.  The  initial  conditions  are 
characterized  by  flight  at  zero  or  small  bank  angles.  During 
the  maneuver,  large  bank  angle  changes  may  be  expected,  together 
vltu  airspeed  changes  that  tend  to  increase  as  the  length  of 
the  maneuver  is  prolonged. 

Accordingly,  the  initial  conditions  selected  for  this  type 
of  study  are  steady  straight  s.'snetric  flight  (see  Definitions 
and  Section  5-2-2),  where  the  following  quantities  are  eero: 

Pi  P1  ^1  R1  *1 

In  the  expansions  used  in  these  equations  of  motion  the  small- 
disturbance  linear  velocity  variables  are  used  and  the  aerodynamic 
forces  and  moments  are  linearized  as  in  Equation  (30).  It  may  be 
appreciated  that  in  particular  cases,  non-linear  aerodynamic 
force  and  moment  ’characteristics  may  be  inserted  into  these 
equations . 

Without  repeating  the  lengthy  algebraic  manipulations  involved, 
the  development  of  the  tracking  equations  of  motion  may  be  out- 
lined in  these  steps: 

1.  Subs  dtute  *he  velocity,  angular  position  and  aerodynamic 
force  ana  moment  expansion  equations  (1C)  to  (15),  (27), 

(20),  and  (30)  into  Equations  (1A)  and  (2A). 

2.  Cancel  the  initial  aerodynamic  force  and  moment  coefficients, 
using  Equations  (6)  and  (7). 

The  final  equations  of  motion  arc  shown  r i matrix  form  in  Ihble  I. 
The  column  headings  of  Table  I are  the  variables  of  motion  and 
Aujctions  such  as  the  products  cf  the  variables  of  nation.  Toe 
rows  represent  separate,  independent  equations  of  motion.  The 
tabular  entries  ore  coefficients  of  the  variables  or  functions  of 
the  variables  which  appear  in  the  particular  equation.  Thus,  the 
firpt  row  of  Table  I Is  interpreted  as: 

(~  - °*  - “In  0*  Q ♦ ♦ XSn  $e  3 0 

2T  2T  TTf 


P»3*  SO 
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If  principal  axes  are  chosen  as  the  arbitrary  body  axes  of 
lfcbl*  I,  • 0 and  • i|j.  . 

6.2  Equations  Suitable  for  the  Calculation  of  Rapid  Rolling 
Maneuvers 

Modern  high-speed  airplanes  are  capable  of  high  rolling  velocities 
starting  from  either  initially  straight  flight  or  from  pullups 
or  pushdowns.  The  inadequacy  of  the  cl/issical  seal 1 -disturbance 
equations  of  airplane  notion  to  calculate  properly  these  notions 
has  become  veil  known. 

These  significant  portions  of  rapid  rolling  naneuvers  generally 
are  of  short  duration  after  the  initiation  of  the  roll,  and 
airspeed  changes  are  consequently  snail  and  of  negligible 
influence  on  the  notion.  Furthermore,  in  conditions  where  non- 
linear inertial  couplings  play  an  important  part  in  the  notions 
and  the  classical  linearised  equations  of  notion  are  least 
applicable,  aerodynamic  and  gravity  effects  are  correspondingly 
less  significant  and  no  re  susceptible  to  drastic  slmpll  float  ion. 

The  initial  conditions  selected  for  this  type  of  study  are 
steady  pitching  or  straight  symmetric  flight  (see  Definitions 
and  Sections  5*2.1  and  5*2.2),  where  the  follovlng  quantities 
are  zero:  Rj.,  P].,  Ri,  . The  snail-disturbance  linear  velocity 

variables  are  used  and  the  aerodynamic  forces  and  moment 9 are 
linearized  as  In  Equations  (30).  Also,  constant  longitudinal, 
velocity  (u  • u » 0),  and  zero  pitch  ang*~,  as  it  affects  the 
force  due  to  gravity  (©i  » • » 0)  are  assumed. 

The  final  equations  of  notion  for  this  case  are  derived  very  ouch 
as  outlined  in  the  previous  section,  and  these  equations  are 
shown  in  matrix  form  in  Thble  II.  Principal  axes  are  chosen 
as  the  arbitrary  body  axes,  in  order  to  reduce  the  required 
multiplications  of  variables.  This  appears  to  be  a reasonable 
simplification  for  practical  computations,  since  there  is  no 
tracking  control  data  system  with  arbitrary  body  axis  orientations 
to  be  supplied  with  computed  airplane  motions,  as  in  the  first 
example. 


t / 


7.0  STABILITY  DERIVATIVES  ON  BODY  AXES 

The  body  axes  stability  derivatives  appearing  in  Tables  1 and  II 
are  related  to  wind-tunnel  data  referred  to  stability  axes  in  this 
section.  It  is  assumed  that  moment  transfers  to  the  assumed  airplane 
center  of  gravity  have  already  been  performed  on  the  wind-tunnel  data# 
if  needed. 


7.1  Longitudinal  Velocity  Derivatives 

The  steps  involved  in  deriving  the  longitudinal  velocity 
derivatives  on  body  axes  fro  a.  wind-tunnel  data  are  somewhat 
tedious#  ond  they  will  be  fully  illustrated  in  only  one  case. 

The  illustrative  example  is  the  body  axes  derivative  xuit  defined 
••  DX  1 . The  X force  may  be  written  as: 

Du'  qxS 

X - C*  (/>/2)V2  S (31) 


Di  f ferentiat  J.ng 

% ■ f V!2  . CXl  (f,  S *1 

Row,  u'  « and  DV/DU  * cos  ca 
Then: 

|3'  s c*i  C0B 


dv 

ffC 

from  page  146  of  Reference  4. 

(32) 


The  derivative  Cx„<  can  have  contrlbutic ns  due  to  compressibility 
ana  aeroelastlc  effects  in  addition  to  the  effects  of  angle  of 
attack  treated  in  Reference  4.  Thus: 


CXU»  • cXa 


Da 

Du* 


DM 

du' 


4 Cy 


d3 

«u’ 


From  page  146  of  reference  4: 


(33) 


Dor  sin  a, 

DU  TI~ 

and 

- - sin  crA  (34) 

Thu  partial  derivatives  DM/0u’  and  Dq/Du'  for  arbitrary 
body  axee  ere  not  treated  in  Reference  4,  so  these  are  derived 
in  the  following  steps: 

m . \ • Vi  * w 

® a 

- Ml  (1  *^[) 


(35) 
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lb*  differential  <f*  la  obtained  from  the  total  velocity  equation: 
• U*  ♦ V2  ♦ W2 
Differentiating: 

2*1  <W  • 2UidU  4 2VxdV  <•>  2WidW 

With  the  initial  lateral  velocity  Vi  equal  to  tero  and  vlth  u'» 
dW/Vl  » Ui/tl  » cos  ax  , and  W/^x  * Bin  c^: 

oq  ♦ a'  sin  oq  (36) 

Equation  (36)  into  Equation  (35),  and  taking 


dU/Vl  , O'  - 

^ ■ u'  cos 

*1 

Substituting 

partlals: 


|2,  . Mx  cos  ox 

IS*  ■ mi  °in  °i 

Again: 

V ■ Vi  (1  .3) 


Pron  Equation  (36),  using  q ■ ( Pr*)\ 2,  and  neglecting  products 
of  u'  and  o'  : 

* • f V (l  4-  2u'  COS  Ox  + 2 a'  #in  °l)  (38) 


Takiug  partial  derivatives: 

lu'  ' ^l2  008  al 

la'  6ln  °1 

In  the  next  part  of  this  derivation,  the  longitudinal  and  normal 
aerodynamic  force  coefficients  C*  and  are  related  to  the 
force  coefficients  measured  on  wind-tunnel  stability  axes  Cp, 

Op,  and  to  the  calculated  net  thrust  coefficient  Tc'.  It  is 
assumed  that  the  effects  of  running  jet  engines  are  all  represented 
with  sufficient  accuracy  in  the  wind-tunnel  tests  except  for  the 
net  thrust.  That  i 0,  it  is  assumed  thuc.  the  intake  normal  force 
discu8oed  in  Reference  5 and  any  air  intake  drag  forces  are 
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included  in  the  forces  and  moments  measured  in  the  wind-tunnel. 
This  vlll  be  the  case  when  the  full-scale  Jet  inlet  velocity 
ratio  is  matched  in  the  wind  tunnel,  and  the  duct  inlet  details 
are  represented.  With  these  conditions,  and  are  obtained 
from  Figure  5 as: 


• Cj^  sin  a - Cp  cos  a + Tc’  cos  ip 

C j • -C^  cos  a - Cp  sin  a - Tc*  sinl^ 

Equations  (40)  can  be  differentiated,  yielding  for  the 
derivatives  of  Equation  (33 ); 


u, 

X - partial 


% “ (CLq  + W ®in  °i  + (Cl1  * C°a)  008  °i 

C«  “ Chi  8in  °i  ■ °I>M  C03  01  + V„  008  1t  <41) 

cXq  * CL.  8in  «L  - °Dq  cos  °i 


The  variation  of  thrust  with  airspeed  represented  by  the  partial 
derivative  is  more  conveniently  expressed  in  terms  of  the 

rate  of  change  of  net  thrust  with  airspeed  Tly . This  is  accomplished 
in  the  following  steps.  Let. 

T * Tc*  (/>/*)  V2 S 

Then,  as  in  the  derivation  of  Equation  (32), 


3T 

Zn' 

' T«'u-  * 

cos  cq 

(42) 

Now: 

3T  . 

0% 

dT  1 ZU 
3u'  V dV 

(43) 

and 

3U/3V  = l/ cos 

cq,  from  Reference  4.  Also: 

Vu 

, = a?c’  dM 
<5  M ou* 

* V,.  M1  C03 

(44) 

Substituting  Equations  (43)  arid  (44)  into  Equation  (42), 

vM-\  h ' f*i  (45) 

Ml 

The  final  step  required  is  the  substitution  of  Equations  ( 33 ) » 
(34),  (37)*  (39),  (40),  (4l),  and  (4$)  into  Equation  (32),  and 
simplifying . When  this  procedure  has  been  carried  out,  the 
results  roy  be  summarized  as  on  the  first  row  of  Tuble  III.  The 
oth'*r  longitudinal  velocity  derivatives  on  body  axes  are  obtained 
in  like  dinner. 
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7.2  Angle  of  Attack  Derivatives 

The  derivation  of  the  angle  of  attack  derivatives  on  body  axes 
is  illustrated  in  a single  ease,  as  for  the  longitudinal  velocity 
derivatives.  The  illustrative  example  is  xa<  , defined  as 
• X 1 . Differentiating  Equation  (31)  with  respect  to  V,  and 

OB’  qjS  w 

noting  that  o'  • ^ and  0V/&W  * ein  0^  from  page  146  of  Refer- 
ence 4, 


dX  1 

da'  q^§ 


cXa.  + 2cx1 


Hov: 

<v 


da 

po' 


+ CY  .♦  C 
XM  aa' 


(**6) 

(•*7) 


From  page  146  of  Reference  4; 

9a  e cos 
»y  — 


and 


Jf.  * COfl  °i  (W) 

Finally,  by  substituting  Equations  (37),  (39),  (40),  (4\),  (45), 
(47),  and  (48)  into  Equation  (46),  the  derivative  xQi  is  obtained 
as  the  first  row  of  Table  IV.  The  reTainlng  angle  of  attack 
derivatives  on  body  axes  are  obtained  in  like  winner. 

7.3  Sideslip  Derivatives 


The  relationships  between  the  sideslip  derivatives  on  arbitrary 
body  axes  and  the  sideslip  ierivatives  on  stability  axe6  obtained 
ir  the  wind  tunnel  are  derived  on  page  61  of  Reference  4.  These 
equations  are  reproduced  in  Table  V of  this  report,  using  modern 
notation.  Second-order  derivatives  like  y^i^,  are  included  in 

the  equations  of  Tables  I and  II,  having  the  significance  of  the 
rate  of  change  with  angle  of  attack  of  the  sides) ip  derivative 
in  question.  Formulae  lor  obtaining  the  second-order  sideslip 
derivatives  directly  from  wind-tunnel  data  are  included  in  Table 
V. 


D 
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Sysatrlc  Rotary  Derivatives 


Ilw  derivation  of  the  symmetric  rotary  derivatives  on  body  axes 
in  terms  of  vind-tunnel  or  calculated  data  on  stability  axes  is 
illustrated  by  the  derivation  of  x£  , defined  as  dX  1 . 

i4  qjfi 

Differentiating  Equation  (3l)  with  respect  to  q and  noting  that 

8V/9%  • 0, 


DX  _1_ 

«q  qxs 


(*9) 


The  disturbance  pitching  velocity  variable  q is  Independent 
of  angle  of  attack,  Mach  number,  dynamic  pressure,  and  the  other 
variables  of  the  problem.  Thus,  differentiating  Equation  (4o) 
with  respect  to  q, 


cXe  * CI^  * x - CDq  co*  (50) 

The  other  symetric  rotary  derivatives  s'uoarlzed  in  Table  71 
are  derived  in  like  manner. 


7.5 


Asynwetric  Rotary  Derivatives 


The  relationships  between  the  rolling  and  ya»  'g  or  asymnetrlc 
rotary  derivatives  on  arbitrary  body  axes  ar . these  same 
derivatives  referred  to  stability  axes  are  . erived  cn  pages  71 
and  72  of  Reference  4.  These  expressions  are  complicated  by 
the  need  to  resolve  not  only  the  forces  and  moments  ou  the 
body  axes,  but  also  the  rolling  and  yavlng  velocities.  The 
results  of  these  derivations  ere  summarized  in  Table  VII, 
including  the  second-order  asymmetric  rotary  derivatives. 


7.6  Control  Surface  Derivatives 


The  control  surface  derivatives  n body  axes  in  terms  of  vind- 
tunnel  data  on  stability  axes  may  be  derived  in  a similar  manner 
to  the  symmetric  rotary  derivatives  in  Section  7.4.  This  is  the 
•nsequence  of  the  independence  of  the  control  surface  distur- 
bance angles  8 with  respect  to  the  other  variables  of  the  motion. 
The  body  axes  control  surface  derivatives  are  shown  in  Table  VIII. 
Second-order  control  surface  derivatives  like  x*.  a d /dX»  1 

*a  da’  <S6>^S 

ore  not  ohown  in  Tuble  I,  II,  or  VIII,  but  these  can  be  readily 
deluded  if  needed. 


V 
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where  [oj  le  a ntrli  representing  the  successive  non- inter- 
acting  rotations  perfbmed  by  the  ease  of  the  neaeuring 
instrument  to  arrive  at  the  airplane '■  orientation  and  £l]  le 
the  orientation  satrlx  of  Reference  3.  7T  the  attitude  -neaeuring 
lustrums nts  are  conventional  vertical  and  directional  free  gyros, 
Squstlon  ($3)  reduces  to  Equations  (lU)  and  (16)  of  Reference  3. 


8.3  Readings  of  Velocity  - Measuring  Instruments 

In  the  most  general  case,  velocity -aeasurlng  instruments  My 
have  any  arbitrary  location  and  orientation  vlth  relation  to  a 
aaneaverlng  airplane.  Of  course,  unless  the  velocity-measuring 
Instruments  are  located  at  the  airplane's  center  of  gravity  and 
are  aligned  vlth  the  arbitrary  ? vs ten  of  boiy  axes,  the  Instrument 
feedings  (corrected  for  local-flow  distortion  position  errors 
such  as  upwash  or  sldeuash,  instrument  errors  and  reduced  to 
true  speed)  will  not  be  given  by  V,  V,  and  W of  Equations  (5), 

(5 A),  or  (19). 


A general  set  of  instrument  axes  nay  be  selected,  vi*\  origin 
at  a radiuc  vector’  r from  the  airplane's  center  o'  gravity, 
and  vlth  arbitrary  orientation  angles  tj,  9^,  » ♦.  vlth  respect 

to  the  body  axis  system.  These  general  4 atruaent  axes  are 
illustrated  in  Figure  6.  The  true  spee/  on  instrument  axes 
nay  be  expressed  as: 


Yx  - Y * 4 l (5*0 

where 


Yj  • L ui  ♦ * Vj  ♦ n Wj 

\ - lU+JV+kV 

A • IF+jQ+kF 

r • Ix+^y+fc* 

Equation  (5*0  may  be  expressed  In  Cartesian  form  using  the 
Orientation  Matrix  [L]  of  Reference  3#  but  with  j/ j,  6j,  and 

In  place  of  ♦ , 9,  and  ♦.  'ftius: 


UI 

U - Ry  ♦ Qz 

V. 

i 

‘ M 

V - Pz  ♦ Rx 

y 

WI 

m m 

V - Qx  ♦ Py 

» a 

(55) 
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Uj  ■ (U  - Fy  ♦ Qx)  cos  0j  cos  tj 

♦ (V  - Pz  * Fx)  cos  ©j  sin  tj 

- (V  - Qx  ♦ Py)  sin  ©j 

Vj  - (U  - Fy  + Q r)  (cob  VjSlnBjslndj  - slnfjcostj) 

♦ (V  - Px  ♦ Fx)  (BLnfx8^n*l8iodI  4 coafjcooBj) 

+ (W  - Qx  ♦ Py)  ala  Bj  cob  ©j 

Wj  • (U  - Fy  ♦ Qx)  (cosfjcosBjstn© j ♦ tlnBjslnBj) 

♦ (V  - P*  ♦ Fx)  (slatfCoesjslndj  * cost^sint^) 

♦ (W  - Qx  ♦ Py)  COB  ©j  COB  Bj 

In  ths  usual  flight  instrument  arrangement  where  sideslip  and 
angle  of  attack  vanes  are  provided  which  are  Independent  of 
cz*  another,  only  the  angle  of  attack  signal  Is  consistent 
vltta  the  conventions  of  Figure  3*  For  this  arrangement,  the 
Instrument  readings  are  (neglecting  local  flow  distortions 
such  as  upwasb  or  sidewash): 

Pj  - tan'1  (Vj/Uj)  J 

-1  | (56) 

Oj  - tan  (Wj/Uj)  I 

Although  Equations  (55)  and  (56)  would  be  cunlersome  for 
computational  purposes,  considerable  simplification  is  usually 
possible  through  tuc  neglect,  v?  «msll  terms.  For  example, 
the  orientation  angles  Bj  and  will  clmost  invariably  be 
equal  to  zero,  one  or  more  of  tfie  distant  *0  x,  y,  and  ?.  will 
usually  be  negligible,  and  the  use  of  the  small-disturbance 
linear  velocity  variables  u',  ©',  and  a 'will  probably  Justify 
eliminating  he  tan"-1  from  Equations  (56). 

8.k  Readings  of  Acceleration-Measuring  Instrument? 


As  in  ‘he  case  of  velocity-measuring  Instruments,  ecccleratlon- 
measurii.g  instruments  may  have  arbitrary  locations  and  orienta- 
tions vjth  respect  to  the  airplane.  A voctor  equation  for  the 
linear  accelerations  on  instrument  axes  may  be  written  as: 

”V_  “ \ + J\  * \ + A-  x r + J\.  x (fi.  x r)  (57) 
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vtett,  In  addition  to  tbe  previous  definitions: 

I,  ■ i #,  * ! *1  * S 4I 

A » tp 

A comparison  of  Equation  ($7)  vitb  tbs  vector  font  of  tbe 
equations  of  motion  of  tbe  accelerometer  and  of  tbe  airplane 
itself  leads  to  an  Important  principle  applying  to  acceleration* 
measuring  instruments.  This  principle  may  be  stated  as 
follows : 


"Tbe  reacts*  along  any  direction  of  a linear 
accelerometer  mounted  at  the  airplane's  center 
of  gravity  is  proportional  to  tbe  resultant 
aerodynamic  force  (including  thrust ) in  that 
direction,  and  is  independent  of  the  airplane's 
attitude  In  space." 

Because  of  tbe  importance  of  this  principle,  It  Is 
derived  from  basic  considerations  in  appendix  A. 


Although  lengthy,  tbe  general  expressions  for  linear  accelerations 
along  Instrument  axes  are  presented  below  in  Cartesian  form. 

These  equations  are  derived  from  tbe  orientation  matrix  ("Lj] 
used  in  Equation  (55)*  Thus:  L J 

x (R2+$2)+y(PQ-R)+z(RP*Q) 

(58) 


► « 

Vi 

■ N 

*i 

4 

or 


* 

U, 


■ [u-RV+CW-x  (R2*^2 )*y ( PQ-R )♦* (RP+Q )]  coaQjCosij 

♦ [v-PW+RU+x  (PQ*R)-y(P2*R2)+z«}R-P)]  cose^in^ 

-[w-QU*rV+x  (PR-Q)*y(QR*P)-s(Q2*P2)]  ainBj 

« J^O-RV+QW-x  (R2*Q2)*y(PQ-R)+z(RP*Q)]  (cosfjSlnSjelnBj-ain^jcosSj) 
+ [v»PW4RU+x  (PQ*&)-y(P2*R2)+z(3R-P)]  (sinfjainSjainQj+cosijCoasj) 

♦ [w-3U+FV«-x  (PR-$)+y(QR*P)-i(Q2+F2)]  fcln»j  coaSj 


•Neglecting  instrument  lag  and  errors. 


Vj  » [u-RV-KH-x  (R®+Q®)+y(PQ^)4*(RP+4)]  (costjcos*jsln6j+sintj8ln# j) 
+[V'PW+RU4x  (PQ-*£)-y(p2-*R2)-«(<|Up)]  (slntgcoss^lnep-costjsln#^) 
♦[w-QU+FV+x  (H»-$)+y('*+£ )-*(Q2+P2)]  co»6j  cestj 

As  in  tbs  ease  of  the  linear  velocities  of  Equations  (55)/  the 
orientation  angles  tj  end  ♦ j will  almost  invariably  be  equal  to 

sero,  and  one  or  store  of  the  distances  x,  y,  and  s (referred  to 
airplane  body  axes)  will  usually  be  either  sero  or  negligibly 

snail. 


Equations  (57)  and  (56)  give  the  linear  accelerations  on  instrument 
axes,  in  vector  and  Cartesian  fora,  respective!,-.  To  go  from  these 
equations  to  the  readings  of  actual  acceleration-measuring  Instru- 
ments requires  consideration  of  three  additional  factors: 

These  are: 

1.  The  effects  of  gravity  forces  on  the  suspended  or  pivoted 
suss  (sometimes  called  the  seismic  element). 

2.  The  dynamics  of  the  Instrument  Itself,  which  is  assumed  to 
have  spring  and  damping  forces  acting  on  the  suspended  mass. 

3.  The  effects  of  angular  accelerations  of  the  Instrument 
mounting,  when  the  accelerometer  pass  1c  pivoted  as  a 
pendulum. 

The  treatment  used  to  derive  the  readings  of  actual  acceleration- 
measuring instruments  follows  generally  that  of  Reference  6.  The 
effects  of  gravity  forces  on  the  suspended  or  pivoted  nase  are 
based  on  the  following  gravity  resolution  along  lnstrua<ont  axes 
(using  12quation(5)of  Reference  3): 


[*i]  • M 


j-g  sin  e 
g sin  ♦ cos  © 
g ccs  ♦ cos  6 

Expanding: 

gXj  • gf-coe*!  co.*x  slno  <•  cos9j  sin^j  sin#  cos  0 
-sinQj  cos*  cos©} 

g^  ■ g £-(sin#j  sinSj  cos|j  - sin*j  coa#x)  sin  © 

♦(sin^j  sinQj  ein*x  + coo#x  coo#j)  sin#  cos© 
+ sin«j  co90j  cos#  cpso] 


(59) 
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ggj  » c [-  (eoBfjOoMjtlnSi  4 sintjslnv  j)  sin  e 

4 (slntjcossjslnej  - ccsfjslnSj)  slo«  cote 
4 cot* jeo*6jooa* cote j 


In  the  application  of  Equations  (56)  and  (59 )#  two  types  of 
linear  acceleration-measuring  lnstnnents  are  distinguished, 
the  linear-type  aeoeleroaeter  and  the  pendulun-type  accelero- 
Mttr  with  a rotary  transducer.  The  sensitive  direction  of  the 
linear-type  accelerometer  la  assumed  to  be  parallel  to  one  of 
the  Instrument  Axes  Xj,  Yj,  or  Zj  of  Figure  6.  The  sensitive 

direction  of  the  pendulun-type  accelerometer  la  defined  by  the 
orientation  of  the  support  axis  parallel  to  one  of  the  Instrvuent 
Axes.  For  sero  signal  output,  the  pendulum  arm,  the  support  axis, 
and  the  sensitive  direction  of  the  accelerometer  form  an  orthogonal 
system.  By  definition,  we  define  the  Instrument  Axes  for  this  case 
as  being  pazmllel  to  each  of  the  orthogonal  directions  formed  by 
the  pendulw  arm,  the  support  axis,  and  the  sensitive  direction. 

The  orientations  of  llnesr-and  pendulum-type  accelerometers 
relative  to  Instrument  Axes  are  Illustrated  In  Figure  7* 


8.4.1  Linear-Type  Accelerometer 


The  transfer  function  of  s linear-type  accelerometer  is  given  by 
(in  the  notation  of  Reference  6): 


where: 


(60) 


T • Accelerometer  output,  or  the  disturbed  motion  of  the 

seismic  element  relative  to  the  case,  parallel  to  either 
the  Xj,  Yj  or  Zj  axes  (figure  7),  t has  the  dimensions 

of  Length. 

s^a#  - Resultant  linear  acceleration  of  the  casf  in  the  sensitive 
direction,  s^s.  Is  equal  to  Oj,  Vj,  or  Wj  of  liquation  (56), 

depending  upon  the  sensitive  direction  of  the  instrument. 
*2aa  lv  parallel  to  the  direction  of  T . 

P « Tilt  ang.le  of  the  sensitive  direction  of  the  accelerometer 
to  the  horizontal,  p is  equal  to  g^j/g,  nr  gfcj/ g 

of  Equstio"  (59)#  depending  upon  the  sensitive  direction  cf 
the  instrument,  and  is  normal  to  f and  s*nfl. 
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J - Dynamic  properties  of  the  linear  accelerometer. 

K - gm'A 
Wn  - Vk/n* 

J - d/2  V»'k 

vbere  «’  is  the  seismic  bass 

k is  the  spring  constant 
d is  the  damp log  constant 

Equation  (60)  assumes  that  the  accelerometer  is  essentially  a 
linear  device.  In  the  sense  that  the  output  T It  directly  pro- 
portional to  the  i.iout  *?ag  ^ at  any  particul  \r  frequency. 

this  assumption  Is  not  at  all  inconsistent  vith  assumptions  of 
non-linearity  for  the  motions  of  the  airplane,  since  It  is  not 
unusual  to  find  accelerometers  with  large  usable  linear  ranges. 

8.4.2  Pendulum-Type  Accelerometer 

A generalized  transfer  function  for  a pendulum- type  accelerometer 
may  be  derived  using  the  methods  of  deference  6.  Three  directions 
are  defined  as  In  Figure  7,  vlth  unit  vectors: 

s • Sensitive  direction 

t ■ Normal  to  sensitive  directlcn  and  support  axis 

u ■ Support  axis 

From  Equation  (4)  of  Keference  6: 

7 £ ■ S - I * a (61) 

where i 

J ■ 0 » + 0 t +Tu 

% 

2 • 0s+0t  + Qu 

ytj  - Ag 8 ♦AyU 

f • *Tfs  + ft  + 0u 
2*  * \ a - Vt  t ♦ u 


I 


K 

» 

f 


I 

t 
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Tbs  moment  applied  to  the  pendu.lv  relative  to  its  case  Q it: 

Q > -dT  - kT  ♦ a f x g (62) 

where: 

i - + 

Substituting  Equation  (62)  into  Equation  (6l)  and  taking  the 
resultant  along  u: 

* . * 1 <«3> 


vbere: 


K 


1 


♦ «t  - it 

6 S 


♦ »,fst 


1 


2 ^/l(k  ♦ m'  fgt  • m' 

The  following  definitions  apply  to  Equation  (63): 

T ■ Accelerometer  - cput  or  disturbed  rotary  motion  of  the 

pendulum  arm  aoout  its  support  a- is  relative  to  the  case. 

T 4a  measured  about  the  support  axis  u , which  is  parallel 
to  one  of  the  Instrument  Axes  Xj, 

V#  » Beaultant  linear  acceleration  of  the  case  in  its  sensitive 
direction,  along  s.  Va  Is  equal  t->  either  Up  7p  ov  !>j 
of  Equation  (58)  depending  upon  the  orientation  of  the 
lnetrument. 

/I  ■ Angular  velocity  of  accelerometer  case  about  the  pendulum 

u cupport  axis  direction  u.  nu  i8  aqual  to  cither  Pp  Qp 
*1* 

ga  ■ Component  of  the  acceleration  of  entity  along  the  acceler- 
ometer sensitive  axle  s.  g6  is  equal  to  either  g^,  gy^,  or 

of  Equation  (5>),  depending  on  the  instrument  orientation. 


- Coeponent  of  tbs  acceleration  of  gravity  along  the  normal 
to  the  sensitive  direction  and  support  axis,  or  along  t. 

• Resultant  linear  acceleration  of  the  case  along  the  norrsal 
to  the  sensitive  direction  and  support  axis,  or  along  t 

Equation  (63)  dcaonstrates  that  the  response  of  a pendulum- type 
accelerometer  used  to  neasure  linear  acceleration  (along  its 


sensitive  axis)  is  dependent  on  the  linear  acceleration  at  right 
angles  to  the  sensitive  direction.  For  example,  a pendulum- type 
accelerometer  vlthout  a spring  used  as  a weight  -downwards  lateral 
accelerometer  will  have  a static  gain  K value  approximately  equal 
to  1/n,  where  n la  the  normal  load  factor.  In  a 2g  puLlup, 
the  seismic  bass  will  travel  half  as  far  In  response  to  a given 
steady  lateral  acceleration  as  It  would  In  straight  flight.  Also, 
the  natural  i*requencyaJQ  is  seen  to  he  increased  In  a pullup  and 

the  dasqping  ratio  is  decreased. 


8.U.3  Pendulum-] 


Jlar  Accelerometer 


The  transfer  function  of  a pendulum -type  accelerometer  used 
an  angular  accelerometer  la: 


K 1 

3 

■33  i? 


where  the  notation  agrees  with  that  for  the  pendulum- type  linear 
accelerooter  of  Section  &.k.2  except  that: 


< • A 

* 2p7T 
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APPENDIX  A 

THE  READ  IRQ  OF  A LINEAR  ACCELEKOMETEB  AT 
tfafe  AIRfrlAlfe*  6 CENft»  OF  &AVITY 


Tbe  important  principle  that  the  reading*  along  any  direction  of  a 
linear  accelerometer  mounted  at  the  airplane’s  center  of  gravity 
la  proportional  to  tbe  resultant  aerodynamic  force  (including  thrust) 
in  that  direction,  and  is  independent  of  the  airplane's  attitude  in 
space  was  stated  in  Section  8.4  of  this  report.  This  principle  is 
derived  in  this  Appendix  from  basic  considerations.  Letting  primed 
values  refer  to  the  accelerometer’s  seismic  mass,  Wevton'a  law  for 
the  rate  of  change  of  linear  momentum  for  both  the  airplane  and  the 
accelerometer  seismic  mass  is: 

P ♦ sg  • m(V  + Ax  V) 

P'  ♦ m'g  - n'(V'  + A xV 

The  accelerometer  is  assumed  to  be  located  at  the  airplane's  center 
of  gravity  so  that  the  vector  r of  Equation  (57)  equals  zero.  Since 
the  accelerometer's  case  is  carried  with  the  airplane,  A'  *A  and 
\*  ■ Y , If  we  neglect  the  relatively  small  velocities  of  the 
seismic  mass  relative  to  Its  case.  Thus: 


P F' 


The  vector  F*  represents  the  external  forces  applied  to  the  accelero- 
meter seismic  mass.  Since  ve  are  neglecting  instrument  lag  or  dynamic 
effects,  F'  arises  entirely  from  the  centering  spring  on  the  seismic 
mass.  In~the  usual  case  this  spring  is  linear  and  the  accelerometer 
reading  is  proportional  to  F' . The  vector  form  of  Equation  (A-2)  showe 
further  that  the  accelerometer  reading  is  proportional  to  the  total 
aerodynamic  force  F in  the  direction  in  which  the  reading  is  made. 

The  absence  of  the~gravity  vector  from  Equation  ft  -2)  proves  that  the 
instrument  reading  is  ii«iependent  of  the  airplane's  attitude  in  space. 

If  the  accelerometer  is  calibrated  so  that  its  reeling  T is  made 
equal  to  unity  when  F'/m'g  = 1.0,  then  Equation  (A-.-)  shows  that  the 
accelerometer  rt  -ds  dir*vtly  the  airplane's  load  fsetot,  or  aerodynamic 
force -to-velght  ratio  tlong  its  sensitive  direction. 
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* Neglecting  instrument  lags  and  errors. 


iBrnfi 


* /8SUMPTI0HS:  l)  Principal  axes  (l^-O) 

2)  Constant  Longitudinal  ’i>: lccitj  (u'»0) 

3)  Gravity  terms  simplified  (0  » 0) 
k)  ♦ « fp  dt 

5)  Initial  synwtric  flight  (Jj_  Pj  ^ ALL  • 0) 

6)  0.-11  perturbations  in  sideslip  and  angle  of  attach.  ,ia 
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TABLE  III 

LQWOn’UDIWAL  VELOCITY  DERIVATIVES  OH  BODY  AXES 

(CL,  Op  and  C*  are  assisted  to  Include  all 
running-engine  effects  except  for  the  net; 
thrust) 


DERIVATIVE  SYMBOL 


STABILITY  AXES  EQUIVALENT 


*1  1 

du*  qjS 


*u»  co*2ai  + ^‘CLo‘  °Di^  •,ln2al 

♦ ^>a+  m1cIm+  S^Lq)  co»  °i  •*«  C*j,  * r %co*  iTco*°i 


_1_ 
Du*  qxS 


(-^Li*  MlcI^|“2<*lcLq)  co82«i  + (ppa*  CLi)  *^2<\ 
♦(Ci^-CDi-M^C^H-Jiq^Cjj^)  coso^  sina^  - — T^  sin  Lp  cos 


Dl  _1__ 
3u'  qxSb 


(8C|,  + MjC|m  +2^1  q)  cos 


DM  1 
Du'  q^Sc 


<MiSi + 2<ic«q + r ^ ~ - 2rc;  ^2)  c°«  °i 

- C sin  ^ 


SI!  i 

£u'  q^Sb 


(2Cn  ♦ Mj^C_  42ql  C ) cos  0^ 
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1ABI£  IV 

AMPLE  OF  ATTACK  DERIVATIVES  ON  BODY  AXES 

(CL,  Op,  and  Cb  are  ass  used  to  include  all 

running-engine  effects,  except  for  the 
net  thrust) 


DERIVATIVE 

SYMBOL 

STABILITY  AXES  EQUIVALENT 

dx  1 
5a'  qp> 

az  _i_ 
55*  qxS 

V 

V 

^clx*  °Da)  c°9\  + 009  °i  *la  °i 

♦ (2C 1 MiClM+2q1CL^)  sln2^  ♦ cos  ij  sin 

CDi)cos2a1  ♦ (*CLl-  Cp^-MjCj^-aqjCj^)  cos  0^  sin 

♦ (-2^-  MicDM*2qicD4^  #ln2oi  * % #ln  *t  ,ln  °i 
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